Entanglement, Elasticity and Viscous Relaxation of Actin Solutions 
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We have investigated the viscosity and the plateau modulus of actin solutions with a magnetically 
driven rotating disc rheometer. For entangled solutions we observed a scaling of the plateau modulus 
versus concentration with a power of 7/5. The measured terminal relaxation time increases with a 
power 3/2 as a function of polymer length. We interpret the entanglement transition and the scaling 
of the plateau modulus in terms of the tube model for semiflexible polymers. 
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Networks of semiflexible macromolecules are major 
constituents of biological tissue. There is experimental 
evidence [Q-Q that certain aspects of biologically impor- 
tant macromolecules, such as DNA and actin, are well 
described by the minimal theoretical model of a semi- 
flexible macromolecule, also known as the wormlike chain 
model. This model represents the polymer as a smooth 
inextensible contour with an energy cost for bending and 
includes ideal flexible chains as a limiting case. The bend- 
ing modulus of the single molecule can be expected to be 
constitutive also for the collective mechanical properties 
of gels and sufficiently concentrated solutions of semi- 
flexible polymers. (Recently, possible contributions from 
twist have also been discussed H.) However, very little 
is known about how semiflexible polymers build up sta- 
tistical networks, and how the macroscopic stresses and 
strains are mediated to the single molecules in such net- 
works. This is also known as the entanglement problem,. 
In this Letter, we report on experiments performed with 
a magnetically driven rotating disc rheometer, which elu- 
cidate some important aspects of the entanglement prob- 
lem. The systems under scrutiny are in vitro polymer- 
ized actin solutions of various concentrations c and av- 
erage polymer lengths L. Actin ||6|] forms large semiflex- 
ible polymers with a persistence length Ip of about 17 
/im ||7|,^ (comparable to typical filament lengths in our 
experiments), and is the most abundant cytoskeletal el- 
ement in most eucariotic cells. We have analyzed the 
transition from the dilute to the semidilute phase (the 
entanglement transition) as a function of polymer length 
and concentration. The data can be interpreted in terms 
of a virial expansion for effective "tubes" . For entangled 
solutions we observed a scaling of the plateau modulus 
G^ versus actin concentration c. This is compared with 
various theoretical predictions [p|-p^. Lastly, we ana- 
lyzed the dependence of the zero shear rate viscosity on 
polymer length, which exhibits a much weaker length de- 
pendence than one would expect theoretically from work 
by Odijk § andDoi |l|. 

Actin was prepared as previously described pq |, and 
purified in a second step using gel column chromatog- 
raphy (Sephacryl S-300). Monomeric actin (called G- 



actin) was kept in G-buffer, consisting of 2 mM Imidazol 
(pH 7.4), 0.2 mM CaCla, 0.2 mM DTT, 0.5 mM ATP, 
and 0.005 vol% NaNs. Polymerization was initiated by 
adding 1/10 of the sample volume of 10-fold concentrated 
F-buffer containing 20 mM Imidazol (pH 7.4), 2 mM 
CaCla, 1 M KCl, 20 mM MgCls, 2 mM DTT, and 5 mM 
ATP. Gelsolin was prepared from bovine plasma serum 
according to Ref. |0], and stored dissolved in G-buffer 
at 4°C for several weeks. The purity of the proteins was 
checked by SDS-PAGE electrophoresis. After staining 
with coomassie blue |18| only one single band was de- 
tected. The mean length of actin filaments was adjusted 
by adding gelsolin to G-actin before initiating polymer- 
ization. According to results by Janmcy et al. p9[ we 
computed the average actin length from the molar ratio 
TAG of actin to gelsolin as L [/im] — rAG/370. All mea- 
surements were done at room temperature (20 ± 0.1)°C. 
Both oscillatory and creep experiments were performed 
with a magnetically driven rotating disc rheometer, as 
described previously ||l^. Care was taken to keep the 
strain below 1% to probe linear response. For oscilla- 
tory measurements the phase shift between exciting force 
and observed oscillation and the response amplitude were 
recorded. From these two parameters the dynamic stor- 
age and loss modulus (real and imaginary part of the 
stress amplitude divided by the strain amplitude) were 
obtained for frequencies 10/211 = 10^"* to 10^ Hz. The 
creep compliance J{t) was obtained for times t = 10~^ 
to 10^ s by applying a sudden step force to the sample 
and recording its strain, which is proportional to J{t). In 
both cases the apparatus was calibrated with purely vis- 
cous liquids of known viscosities. A quantitative measure 
for the elastic character of a material is the phase shift. 
In the limiting case of a purely elastic medium the phase 
shift is zero, in the opposite case of a purely viscous liquid 
the phase shift is 7r/2. Consequently, the sample behaves 
most rubber-like when the phase shift becomes minimal. 
Therefore, in oscillatory experiments with actin/gelsolin 
the value of the storage modulus at the frequency cor- 
responding to the minimum phase shift was identified as 
the plateau modulus G'^ . For actin samples without gel- 
solin, where no minimum in the phase could be observed 



within the measured frequency range, the storage modu- 
lus at a fixed frequency in the plateau regime was taken 
as G° . This does not affect the functional form of G^ (c) 
but its absolute value. (As a consequence, the vertical 
shift between the two data sets shown in Fig. |^ has no 
physical significance.) For the circles in Fig. 07^° was 
determined by the minimum phase prescription at the 
highest concentration only, whereas relative shifts of G'^ 
at lower concentrations were determined by rescaling to 
superimpose the moduli. The zero shear rate viscosity tjq 
was obtained from measurements of the creep compliance 
J{t) according to rj^^ := linit^oo J{t)/t [0- From creep 
experiments we extracted the frequency dependent mod- 
uU (to obtain G" for Fig. |i|) as described in Ref. ^. It 
was checked that the results agree well with correspond- 
ing oscillatory measurements ||21| . 

Figs. |l|, H, and g show the plateau modulus G° as a 
function of filament length L and actin concentration c, 
respectively. The data in Fig. ^ clearly indicate a transi- 
tion with increasing length of polymers. Similar results 
also have been obtained by Janmey et. al. |22| ] , recently. 
At first sight, one might be tempted to attribute this 
transition to the mutual steric hindrance in a solution of 
rods at the overlap concentration. The observed transi- 
tion is indeed in a parameter regime, where the polymer 
length L is not much larger (about a factor of 5) than the 
mesh size ^m, and we originally attempted to interpret 
the data this way. Some more thought suggests, however, 
that there is no transition expected for the plateau mod- 
ulus of stiff rods; a sudden increase in the shear modulus 
near the overlap concentration would not be in accord 
with the virial expansion for the osmotic pressure of rods 
p3[ , which predicts a smooth dependence on c and L 
below the nematic transition. One can hardly imagine 
the shear modulus of a semidilute solution of rods to be 
larger than the osmotic compression modulus. The so- 
lution could easily escape the shear stress by local com- 
pression. On the other hand, the actin solutions in our 
experiments were below the critical concentration for the 
nematic transition |^,^ . We can thus conclude that the 
observed sudden increase (Fig. |lh , and the enhanced con- 
centration dependence (Fig. g,^ of G° above a certain 
threshold are related to the semifiexible nature of actin 
filaments. Their persistence length of about 17 /^m [0,H, 
albeit much larger than the typical mesh size ^m — 0.1 — 1 
/im of the studied networks, can not be assumed to be in- 
finitely large. Otherwise the data would have to obey the 
prediction of the classical theory for dilute and semidilute 
rods [^ 



G" = Si/fcsT/S 



(1) 



Here v = 3/£,'^L is the polymer number density. It is not 
conceivable that the sudden steep increase of G" with 
polymer length is merely due to internal modes of single 
polymers neglected in the theory for stiff rods. Instead 
we are forced to look for a cooperative effect. 
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FIG. 1. The plateau modulus above the entanglement tran- 
sition as a function of polymer length for constant monomeric 
actin concentration c — 1.0 mg/ml. The solid line corresponds 
to Eq. (H) with ^i =0.38. The increase of G° for large L is 
not yet fully understood. 
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FIG. 2. The plateau modulus near the entanglement transi- 
tion as a function of polymer concentration for short rod-like 
actin filaments {L = 1.5 /im). Two different methods were 
used to extract G" from the data (see main text). Also shown 
is the theoretical prediction, Eqs. (bl) for ^i — 0.47. Theoret- 
ically the transition is expected at c* — 0.68 mg/ml. 
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FIG. 3. Concentration dependence of the plateau modu- 
lus of pure actin (□) and actin with a small amount of gel- 
solin (rAG ~ 6000 : 1) corresponding to an average actin 
filament length L = 16 ^m (O). The straight lines indicate 
the power 7/5 corresponding to the scafing limit of Eq. (H) 
with ^1 — 0.40 and ^i = 0.46, respectively. 



In the following we attempt to give a simple interpre- 
tation of our observations based on the tube concept for 



semiflexible polymers developed theoretically by Odijk 
|g[ and Semenov [Q and related to the shear modu- 
lus by Isambert and Maggs |12], recently. Experimen- 
tally, it has been demonstrated by video-microscopy that 
in semidilute actin solutions the filaments are confined 
to tube- like cages [^. These cages severely hinder not 
only transverse and rotational motions but also undu- 
lations on length scales larger than a certain length Lg, 
called deflection length or entanglement length. Using the 
wormlike chain free energy one can relate Lg to the tube 
diameter d hj L^ = \/2d^£p P,p8[|. On the time scale 
of the plateau, modes of wavelength smaller than Lg are 
already equilibrated. Hence, on this coarse grained time 
scale we can think of the polymer solution in terms of 
an ensemble of tubes. If we apply a reasoning similar to 
that used for the osmotic pressure of dispersed rods ||2^ ] 
to the plateau modulus of dispersed tubes of length L and 
diameter d, we can replace Eq. (|l|) by a virial expansion 



ues Le « 0.58^1^^4^^ and d « 0.37^. 
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G° ^ {3/5)ukBT{l + 2B2i'...) 



(2) 



Here, B2 is the second virial coefficient for the tubes. 
Higher order terms are negligible for small volume frac- 
tion of the tubes. (The latter turns out to be less than 
0.1 in our case.) The product 321^ counts the average 
number of collisions of the tubes, and can thus be used 
to define a collision length Lc '■— L/2B2V (always two 
tubes are involved in a collision). According to Ref. |23] 
the second virial coefficient is given by the excluded vol- 
ume B2 = TTdL'^/4 of a tube. However, to stay consistent 
with our assumption that short wavelength modes have 
already relaxed, we subtract from L half the collision 
length at each end to account for the reduced efficiency 
of dangling ends in the entanglement process. The above 
relation between the second virial coefficient and the col- 
lision length thus becomes 



L/Lc - 1 = nvdl^L - Lcf/2 



(3) 



We can determine the still unknown tube diameter d from 
the following consistency requirement. Following On- 
sager's argument for the second virial coefficient we have 
to pay a price in free energy of the order of fc^T per length 
Lc to add a new tube to the solution. On the other hand, 
to suppress thermal undulations of wavelengths larger 
than Lg the tube has to supply a confinement energy of 
the order ksTL/Li. to the enclosed polymer. Now, if we 
want the tube to be a pertinent effective representation of 
the medium surrounding a test polymer in the entangled 
polymer solution, these two energies should be equal. We 
do not actually have to introduce a physical tube into the 
solution when adding a polymer. Hence, for consistency 
we require Lc = L^, i.e., the number of mutual collisions 
of the tubes must equal the number of collisions of the 
polymers with their tubes. For entangled solutions we 
thus find G^ = 9kBT/5^^j^Le, where far from the entan- 
glement transition L^. and d take their asymptotic val- 



The scal- 
ing behavior was predicted by Isambert and Maggs |lj] 
from a different reasoning before. We also note that it 
is included as a limiting case in a more detailed analy- 
sis concerned with the calculation of the absolute value 
of G° [^. The corresponding scaling G^{c) oc c^^^ of 
the plateau modulus is indicated by the solid lines in 
Fig. y. A much stronger concentration dependence - 
as predicted by a purely mechanical model |13| ] or by a 
model with thermodynamic buckling [ pT[ - and the scal- 
ing predicted in [Q are not in accord with our data. On 
the other hand, the agreement of Eq. (||) with the data 
seems to hold beyond the scaling limit of strong entangle- 
ment. To relate the apparent (theoretical) volume frac- 
tion to the nominal experimental actin concentration c we 
introduce the symbol ^1 for the apparent mesh size ^m 
[/im] of a solution with c = 1 mg/ml. Solving Eq. (m we 
predict the entanglement transition to occur at a concen- 
tration c* [mg/ml] = 7 • 2i/4£y2^2/3^(-52^/7)5/2 (.^^gakly 
bending rod limit assumed). The critical concentration 
is thus theoretically by a factor of c*/c « 2.0{ip/L)^/'^ 
larger than the overlap concentration c. For the persis- 
tence length we assumed ip = 17 fj,in [0,|| for all our 
fits. The only free parameter of the theoretical curves in 
Figs.Jl], ||, and ^ is thus ^1 . It was chosen as ^1 = 0.38 for 
Fig. |l|; 6 = 0.47 for Fig. |, and ^ = 0.40/0.46 for the 
upper/lower line in Fig. 0, in reasonable agreement with 
the value ^1 — 0.35 obtained independently by FRAP 
(fluorescence recovery after photo bleaching) pO]. The 
scatter in the value for ^1 merely reflects the poor ex- 
perimental reproducibility of absolute values of G*' for 
F-actin solutions. 

Simultaneously with the length dependence of the 
plateau modulus shown in Fig. |l|, we have measured the 
length dependence of the zero shear rate viscosity 770 . The 
latter is partly due to static effects, namely the length 
dependence of the plateau modulus discussed above, and 
also to dynamics. The terminal relaxation time Tr, the 
characteristic time scale at which a polymer solution be- 
gins to flow, can be obtained up to a numerical coeffi- 
cient from the viscosity via r^ ~ rj^/G^ . Fig. ^ presents 
such data on the length dependence of r^. Data (not 
shown) obtained directly from the frequency dependent 
viscoelastic moduli by the condition G'(27r/r,.) — G°/2 
or by dG" {2tt / Tr) / duj = fall onto the same curve if 
multiplied by numerical prefactors 1.0 and 2.4 [pl| , re- 
spectively. The mechanism for the terminal relaxation 
seems obvious from the tube picture described above. 
Viscous relaxation only occurs when the polymers have 
time to leave their tube-like cages by Brownian motion 
along their axis. The reptation model that was origi- 
nally formulated for flexible polymers, was extended to 
semiflexible chains by Odijk H and Doi |15|. These au- 
thors calculated the disengagement time Td for a semi- 
flexible chain diffusing out of its tube. However, the 



data for r^ presented in Fig. ^ are not in accord with 
their result for t^. The dependence of the observed ter- 
minal relaxation time r^ on polymer length L is substan- 
tially weaker than predicted for r^, even in the stiff limit 
where r^ = lpL/4:D\\ oc L^ (dot-dashed hue in Fig. |^), 
_D|I — ksT /2TTrjL being the longitudinal diffusion coef- 
ficient of the chain in the free draining approximation. 
Instead, the solid line in Fig. H corresponds to the scal- 
ing law Tr oc i"^/^. 
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FIG. 4. Terminal relaxation time in the entangled phase 
for constant actin concentration c = 1.0 mg/ml. The solid 
line indicates the power 3/2 and the dashed line is Eq. (M) 
with a numerical prefactor of 0.10. The dot-dashed line is the 
disentanglement time calculated by Odijk and Doi [|e|,|15{. 

A tentative interpretation of the data can be given in 
terms of a semiflexible polymer diffusing along a strictly 
one dimensional path; i.e., not being allowed to choose 
between infinitely many new directions at its ends. The 
characteristic decay time for self-correlations of the end- 
to-end vector (R(t)R) is then given by p8[ 



TR^ L 



x/D\\{Rr 



{L + 2ip)^/AD\\ 



(4) 



This presents an upper bound for the terminal relaxation 
time within the tube model. As seen from the dashed line 
in Fig. ^, Tfi (for ip = 17 /im) is in fact by a factor of ten 
too large compared to the data but describes fairly well 
the length dependence of r^. The restriction to one path 
implies a very slow decay of conformational correlations 
p8| . An unusually slow decay of stress (the frequency de- 
pendence of G'(w) is still less than linear in the measured 
frequency range) is indeed observed, but this might also 
in part be due to the broad length distribution of actin 
[|l9| . Clearly, further investigations are necessary to come 
to a better understanding of the terminal regime. 

In summary, we were able to measure some important 
physical properties of semiflexible polymer solutions with 
a rotating disc rheometer. We investigated the plateau 
modulus and the zero shear rate viscosity of semidilute 
actin solutions. At a certain concentration c* larger than 
the overlap concentration c an increase in the plateau 
modulus was observed. We interpreted this entanglement 
transition as well as the concentration dependence of the 



plateau modulus in terms of a tube model that takes into 
account the semiflexible nature of the molecules. For 
strongly entangled solutions our data can be character- 
ized by the scaling law G^ ex c^/^ . We also found a power 
law dependence of the terminal relaxation time on poly- 
mer length Tr oc i"^/^, which is substantially weaker than 
predicted for the disengagement time by Odijk and Doi. 
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